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Formaldehyde:  Electronic  structure  calculations  for  the  50 
and  T 1  states8’ 

George  F.  Adams  and  Gary  D.  Bent6*1 

ARRADCOM.  U.S  Army  Ballistic  Research  Laboratory.  Aberdeen  Proving  Ground,  Maryland  21005 

Rodney  J.  Bartlett  and  George  D.  Purvis 

Batclle  Columbus  Laboratory.  Columbus.  Ohio  45201 

(Received  4  June  1980;  accepted  9  April  198 1 ) 

Many-body  perturbation  theory  (MBPT)  and  coupled-cluster  method  calculations  are  reported  for  the 
S„(Y  ‘A  ,|  and  T,(a  ‘A’)  electronic  states  of  formaldehyde.  The  structural  parameters  for  the  Su  minimum 
[Rm  =  1.102  A,  R^o  =  1.211  A,  HCH  =  116.2")  and  the  7",  minimum  (J?tH  =  1.085  A,  Rco  =  1.327  A, 

HCH  =  118",  "out-of-plane”  angle  =  37‘  12')  agree  well  with  experimentally  deduced  values.  Calculated 
heats  of  reaction  for  dissociation  to  radical  products  and  molecular  products  agree  well  with  literature  values. 

The  energy  barriers  for  dissociation  to  molecular  products  and  rearrangement  to  hydroxycarbene  are 
presented.  Vertical  and  adiabatic  transition  energies  are  reported  for  S0  -+  T\ .  while  a  vertical 
transition  energy  for  S.:  *S  .  is  reported. 


I.  INTRODUCTION 

Formaldehyde  is  an  important  and  ubiquitous  chemical 
species.  In  recent  years  the  photochemistry  resulting 
from  excitation  to  the  first  excited  singlet  state  S,  has 
been  the  subject  of  a  number  of  experimental1-4  and 
theoretical5-1  reports  which  raise  a  number  of  impor¬ 
tant  questions  about  interpretation  of  experimental  data. 
Formaldehyde  also  occurs  as  a  crucial  species  in  mod¬ 
els  that  describe  the  flame -zone  chemistry  pertinent  to 
the  combustion  of  carbon-containing  molecules.8  In 
fact,  Fifer  concludes  that  the  vapor -phase  combustion 
of  nitrate  esters  is  dominated  by  the  reaction  sequence 
that  describes  the  oxidation  of  formaldehyde  by  nitro¬ 
gen  dioxide. 8 

Numerous  reports  of  theoretical  studies  on  various 
electronic  states  of  formaldehyde  occur  in  the  litera¬ 
ture.  Recent  papers  by  Goddard  and  Schaefer  (GS),5 
Goddard,  Yamaguchi.  and  Schaefer  (GYS), 7U>  and  Bell8 
provide  excellent  results  as  well  as  thorough  discus¬ 
sions  of  previous  molecular  obital  calculations. 9,10  In 
particular,  the  results  in  the  GS  paper  would  support  a 
tunneling  mechanism11  for  the  decomposition  of  H2CO  to 
H2-  CO.  However,  the  question  of  whether  this  is  re¬ 
quired  to  interpret  the  experimental  results  is  based 
upon  energy  differences  of  only  4  or  5  kcalmol'1.  Since 
the  Cl  calculations  of  GS5  and  GYS1U1  estimate  the  ef¬ 
fects  of  unlinked  diagrams  associated  with  the  Cl  qua¬ 
druple  excitations  by  the  formula  of  Davidson, 12  and 
since  the  estimate  amounts  to  approximately  3  kcai  mol"1 
for  formaldehyde,  it  seems  pertinent  to  study  the  for¬ 
maldehyde  decomposition  with  many-body  methods  that 
properly  treat  the  effects  of  unlinked  diagrams  and 
quadruple  excitations. 13-16  Consequently,  dissociation 
energies  and  energy  barriers  are  predicted  by  MBPT/ 
CCD  for  the  dissociation  to  radical  products 


“'This  research  was  sponsored  in  part  by  the  It.  S.  Army  Re¬ 
search  Office  under  grant  No.  IJAAG29-S0-C-U105. 
#,NAS-NRC  Resident  Research  Associate. 

01  Present  address:  Department  of  Physics,  l  niversity  of 
Connecticutt,  Storrs,  Conn.  iMiZfiH. 


CHjO  —  HCO-r  H  ,  (I) 

molecular  products 

CH20-H2  +  C0  ,  (II) 

and  the  rearrangement  to  trans -hydroxycarbene 

ch2o-hcoh  .  (Ill) 

Other  important  questions  pertain  to  the  electronic 
excitations  in  formaldehyde.  Equations -of-motion  cal¬ 
culations  for  a  number  of  excitation  energies  have  been 
reported. 11,18  We  report  here  MBPT/CCD  calculations 
for  the  ground  (Jf'/l|,  S„)  and  first  excited  {a3A” ,  T,) 
electronic  states  of  formaldehyde.  The  results  include 
equilibrium  structural  parameters  for  each  state,  the 
adiabatic  excitation  energy  for  S„-  T{,  and  vertical 
excitation  energies  for  S0-  T ,  and  S„-S,. 

II.  OUTLINE  OF  THE  CALCULATIONS 

For  the  closed-shell  states  included  in  this  study, 
the  Hartree-Fock  representation  of  a  state  is  given  by 
a  restricted  Hartree-Fock  (RHF)  wave  function,  where¬ 
as  open-shell  states  are  described  by  unrestricted 
Hartree-Fock  (UHF)  wave  functions.  The  integrals 
over  atomic  functions  were  computed  using  Dunning’s 
4s3 p  contraction  of  Huzinaga’s  9s5 p  Gaussian  primitive 
set  for  first-row  atoms,  and  Dunning’s  scaled  (f  =  1.  2) 

3s  contraction  of  Huzinaga’s  4s  primitive  set  for  hydro¬ 
gen.  19,20  A  single  set  of  rf-type  polarization  functions 
augments  the  atomic  basis  sets  for  carbon  and  oxygen, 
with  orbital  exponents  0.75  and  0.85,  respectively.21 
A  set  of  p-type  polarization  functions,  with  a -1.0,  aug¬ 
ments  the  hydrogen  basis  set.21  Integral  calculations 
for  the  MBPT/CCD  calculations  were  performed  using 
the  Moi  H  Ul.i  integral  program22;  electronic  structure 
calculations  were  performed  using  thec.RNiNi  and  pmhi’1 
programs.23 

Structural  parameters  for  the  formaldehyde  con- 
formers  were  obtained  in  several  ways.  Optimal  pa¬ 
rameters  for  the  equilibrium  positions  on  both  the  S0 
and  T,  surfaces  were  determined  by  fitting  the  results " 
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(c)  (d) 

FIG.  1.  Definition  of  structural  parameters  for  formaldehyde 
calculations:  (a)  equilibrium,  (b)  molecular  products  transi¬ 
tion  state,  (c)  rearrangement  transition  state,  (d)  equilibrium 
for  T,  state. 


of  frozen-core,  double -excitation  MBPT  (4)  calcula¬ 
tions.  Such  a  procedure  is  impractical  when  seeking 
transition  state  structures.  Structural  parameters  for 
saddlepoints  on  the  S0  surface  resulted  from  analytic 
gradient  based  calculations.  In  particular,  we  used 
the  gradscf  computer  codes. 24  In  order  to  be  consis¬ 
tent  in  computing  energy  differences  on  the  S0  surface, 
equilibrium  structural  parameters  for  formaldehyde 
were  also  determined  using  the  cradscf  codes.  Finally, 
vibrational  frequencies  were  estimated  using  the  gradi¬ 
ent  codes.  The  same  basis  set  was  used  for  all  calcu¬ 
lations.  The  structural  parameters  for  each  configura- 


TABLE  I.  Structural  parameters  for  formaldehyde  (S„).a 


Equilibrium  configuration 

r 

s 

a 

SCF 

This  work 

1.092 

1. 182 

115.  9 

DZb 

1.084 

1.217 

116.8 

DZ  +  I* 

1.096 

1. 189 

116.3 

DZ  +  P  CIC 

1.100 

1.212 

116.2 

D-MBPT(4) 

1.102 

1.211 

116.2 

Experiment" 

1.099 

1.203 

116.5 

Molecular  products  transition  state 

Ft 

n 

5 

O’ 

t 

SCF 

This  work 

1.606 

1.097 

1.141 

113.3 

1.219 

DZ  +  I* 

1.586 

1. 104 

1.151 

112.5 

1.203 

DZ  *  P  Cl 

1.572 

1.096 

1.179 

110.8 

1.213 

Rearrangement  transition  state 

n 

*-2 

5 

a 

0 

SCF 

This  work 

1.213 

1.095 

1.266 

56.4 

116. G 

DZ  -  P” 

1.262 

1.087 

1.316 

55.0 

117.3 

DZ  •  P  Cl 

1.232 

1.106 

1.303 

55.2 

115.0 

•t'nits:  angstroms,  degrees. 

"Reference  7(a). 

"Reference  5. 

"Reference  24. 

TABLE  11.  Electronic  structure  calculations 
for  CH20(.X  'A,).* 


Level  of  theory 

Structure 

SCF 

D-MBPT<4) 

SCF 

-  113.8990 

-  113. 8974 

D-MBPT(2) 

-114.2672 

-114.2692 

D-MBPT<3) 

-  114.2730 

-114.2738 

D-MBPT(4) 

-114.2818 

-  114.2832 

SDQ-MBPT(4) 

-  114.2814 

-114.2830 

CCD 

-  114.2760 

-  114.2772 

*L'nits:  hartree. 


tion  are  depicted  in  Fig.  1.  Theoretical  and  experi¬ 
mental  structural  parameters  for  each  configuration 
are  listed  in  Table  I.  We  report  dissociation  energy 
and  energy  barrier  results  for  several  levels  of  com¬ 
putation:  SCF;  fourth-order  MBPT  including  all  dou¬ 
ble-excitation  diagrams  D-MBPT(4);  fourth-order 
MBPT  including  all  single-,  double-,  and  quadruple- 
excitation  diagrams  SDQ-MBPT(4);  and  coupled  clus¬ 
ter  double  calculations,  CCD  based  upon  the  wave  func¬ 
tion  exp(r2)<I>0.  Excitation  energies  are  reported  at  both 
SCF  and  MBP^CCD  levels  of  computation.  All  re¬ 
ported  energies  include  the  core  electrons. 

To  compute  dissociation  energies  and  energy  barriers 
requires  accurate  calculation  of  energy  differences.  In 
addition  to  the  electronic  energy,  vibrational,  rotation¬ 
al,  and  translational  energy  differences  must  be  esti¬ 
mated.  We  assume  that  the  rotational  and  translational 
degrees  of  freedom  are  classical.  Vibrational  frequen¬ 
cies  for  the  various  formaldehyde  conformers  were  pre¬ 
dicted  using  SCF-gradient  calculations. 24 

III.  THE  S0  '-4 1 )  STATE  OF  FORMALDEHYDE 

At  its  ground  state  equilibrium  configuration,  formal¬ 
dehyde  possesses  C2„  symmetry.  The  structural  pa¬ 
rameters  obtained  by  several  theoretical  calculations 
are  compared  with  experimentally25  derived  values  in 
Table  I.  The  agreement  between  theory  and  experi¬ 
ment  is  excellent,  especially  for  the  Cl  calculations7<*’ 
and  the  MBPT  results. 

To  determine  chemical  energy  differences  requires  a 
high-quality  ab  initio  calculation  of  the  electronic  ener¬ 
gy  of  those  species  involved  in  the  chemical  reaction. 
The  results  of  several  levels  of  electronic  structure 
calculation  for  the  S0(X  'A,)  state  of  formaldehyde  are 
listed  in  Table  II.  The  Hartree-Fock  energy  for  the 
structure  determined  using  the  gradient -based  SCF  cal¬ 
culations  is  lower  than  that  corresponding  to  the  struc¬ 
ture  found  by  fitting  the  results  of  D-MBPT(4)  calcula¬ 
tions.  However,  through  fourth  order,  the  MBPT  cal¬ 
culations  yield  lower  energies  for  the  structure  obtained 
using  the  D-MBPT  (4)  structure.  In  addition,  the  re¬ 
sults  of  CCD  calculations  predict  a  lower  energy  for 
the  structure  determined  using  the  D-MBPT  (4)  struc¬ 
ture.  Each  of  the  SCF  energies  reported  here  is  a  few 
millihartrees  lower  than  the  single-configuration 
DZ+  P  result  reported  by  Goddard  and  Schaefer.5 
We  expect,  therefore,  that  our  basis  set  is  comparable 
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in  quality  to  that  used  by  GS  in  their  initial  study  of  the 
S„  surface  of  formaldehyde.5  The  two  fourth-order 
MBPT  calculations,  D-MBPT(4)  and  SDQ-MBPT(4), 
agree  almost  exactly  with  one  another.  Clearly  the 
contribution  of  the  single-  and  quadruple-excitation 
diagrams  almost  cancel  at  fourth  order.  The  fourth- 
order  triple  excitation  contributions  would  lower 
the  electronic  energies.  The  CCD  calculation, 
which  includes  higher  order  effects  of  double  and  qua¬ 
druple  excitations,  but  neglects  single-  and  triple-ex¬ 
citation  effects,  predicts  a  smaller  contribution  of  elec¬ 
tron  correlation  than  do  the  fourth-order  MBPT  calcu¬ 
lations.  This  is  due  to  the  neglect  of  single-  and  triple¬ 
excitation  effects,  both  of  which  are  negative.  Although 
CCD  is  not  rigorously  variational,  since  the  principal 
terms  neglected  in  the  calculation  are  negative,  it  is 
quasivariational  and  yields  the  lowest  energy  to  date 
for  the  S0  state  of  formaldehyde. 

IV.  DISSOCIATION  REACTIONS  OF  FORMALDEHYDE 

Thermal  dissociation  of  formaldehyde  is  assumed  to 
produce  radical  products,  hydrogen  atom  and  formyl 
radical.26  Interestingly,  the  dissociation  to  molecular 
products  H2  and  CO  is  approximately  thermoneutral, 
while  the  observed  thermal  dissociation  energy  is  ap¬ 
proximately  88  kcat  mol'1. 21  Both  sets  of  products  cor¬ 
relate  by  symmetry  with  the  ground  electronic  state. 

The  first  singlet  excited  state  of  formaldehyde 
[S,(A  ‘A")|,  has  an  excitation  energy  threshold  lower 
than  the  threshold  for  thermal  dissociation. 1  Photo¬ 
chemical  studies  of  reactions  subsequent  to  excitation 
of  formaldehyde  to  the  state  are  interpreted  in  terms 
of  the  following  three  processes1: 

H2CO(St,  t'')  — H2  +  CO  ,  (IV) 

-H  +  HCO  ,  (V) 

-H2CO(S0,  i/')  +  fev  .  (VI) 

Most  investigators  agree  that  the  formation  of  molecular 
products  (IV)  dominates  for  excitation  energies  near  the 
S,  threshold  (3.495  eV,  80.5  kcalmol'1).  The  impor¬ 
tance  of  the  radical  products  channel  (V)  increases  with 
increasing  excitation  energy.2  For  an  excitation  ener¬ 
gy  equal  to  91. 3  kcal  mol'1,  the  sum  of  quantum  yields 
for  processes  (IV)  and  (V)  is  almost  one,  with  the  molec¬ 
ular  quantum  yield  0. 32 ±  0.03  and  the  radical  quantum 
yield  0.68±0.03.2 

The  crucial  datum  in  understanding  the  photochemical 
dissociation  of  formaldehyde  is  the  energy  difference 
between  the  S,  threshold  (80.  5  kcal  mol'1)  and  the  poten¬ 
tial  energy  barrier  to  the  formation  of  Hj+CO  on  the 
S„  surface.  The  appearance  of  carbon  monoxide  is  much 
slower  than  the  decay  rate  for  S,,  and  the  formation  of 
carbon  monoxide  requires  a  collision.  This  behavior 
suggests  that  the  high-pressure  {>0.1  Torr)  photochem¬ 
istry  of  formaldehyde  excited  to  S,  involves  a  transition 
to  an  intermediate  species  or  to  a  different  electronic 
state  of  formaldehyde.  The  identity  of  the  intermedi¬ 
ate^)  remains  a  point  of  contention,  but  the  results  of 
the  Cl  calculations  by  GS5  and  GYS71*1  imply  that  the  in¬ 
termediate  occurs  on  the  S0  surface. 


Goddard  and  Schaefer  reported  extensive  configura¬ 
tion  interaction  calculations,  including  all  single-  and 
double  excitations  relative  to  their  reference  functions, 
that  elucidated  many  features  of  the  potential  energy 
hypersurface  that  pertain  to  the  photodissociation  of 
formaldehyde. 5  Subsequently,  GYS71*1  reinvestigated 
the  S0  surface  using  analytic  gradient  methods  that  in¬ 
cluded  correlation  effects.  The  calculations  addressed 
reactions, (IV)  and  (V)  and  the  isomerization  of  formal¬ 
dehyde  to  form  trans -hydroxy carbene,  all  on  the  S0  po¬ 
tential  energy  surface.  The  SD-CI  study5  provided  im¬ 
proved  values  for  energy  barriers  to  photochemical  re¬ 
actions  of  formaldehyde.  A  significant  result  of  the 
study  was  the  prediction  that  the  energy  required  for 
isomerization  to  hydroxycarbene  is  only  slightly  greater 
than  the  energies  needed  for  radical  or  molecular  pho¬ 
todissociation.  5  The  SD-CI  calculations  also  predicted 
a  small  energy  barrier  above  the  S,  origin  to  any  pho¬ 
tochemical  reaction  of  formaldehyde. 

In  subsequent  research  Goddard,  Yamaguchi,  and 
Schaefer71*1  reported  features  of  the  formaldehyde  po¬ 
tential  energy  surface  using  analytic  configuration  in¬ 
teraction  techniques.28’29  The  results  of  that  research 
differ  little  from  the  results  obtained  by  GS.5  A  note¬ 
worthy  conclusion  of  the  GYS  work  is  that  the  results 
of  Cl  calculations  using  the  SCF  optimized  geometry 
differs  little  from  those  obtained  using  the  Cl  deter¬ 
mined  transition  state.  For  a  double  zeta  plus  polariza¬ 
tion  basis  set,  the  two  predictions  of  the  molecular 
products  dissociation  barrier  differ  by  only  0. 8 
kcal  mol'1. 


Miller11  employed  the  barrier  height  and  vibrational 
frequencies  predicted  by  GS5  to  calculate  the  collision¬ 
less  decay  rate  for  S0-  H2  to  CO  as  a  function  of  ener¬ 
gy.  The  results  of  Miller’s  analysis  suggest  that  S0 
can  decay  rapidly  to  molecular  products  even  at  ener¬ 
gies  as  much  as  10  kcal  mol'1  lower  than  the  barrier 
height.  Since  that  analysis,  however,  additional  stud¬ 
ies  of  the  S0  surface  by  Harding  et  al. , 50  Goddard  el 
al.  j71*1  and  this  work  have  modified  the  original  esti¬ 
mate  of  the  vibrational  frequencies  for  the  rearrange¬ 
ment  transition  state.5  The  effect  of  these  studies  is 
to  reduce  the  barrier  height  for  the  rearrangement  of 
formaldehyde  to  trans -hydroxycarbene  relative  to  the 
barrier  for  molecular  products  dissociation. 

A.  CH20-CHO+H 

The  theoretical  investigation  of  the  formyl  radical 
using  MBPT  and  CCD  has  been  published.51  We  report 
in  Table  III  the  results  of  those  calculations  that  per¬ 
tain  to  the  study  of  the  dissociation  of  formaldehyde. 
The  experimental  zero-point  vibrational  energies  of 
formaldehyde  and  formyl  radical  are  16.  2”  and  7.8 


TABLE  Ill.  Electronic  structure  results  for  HCOfX  JA').* 


UHF  D-MBPT14)  SDQ-MBPT<4) 

CCD 

Energy 

-113.2769  -113.6343  -113.6364 

-113.6288 

‘Units: 

hartree. 
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TABLE  IV.  Heat  of  reaction  for  the  dissociation  of 
formaldehyde  to  radical  products.* 


Level  of  calculation 

SCF 

77.0 

-8.4 

68.6 

D-MBPT(4) 

94.8 

-8.4 

86.4 

SDQ-MBPT(4) 

93.4 

-8.4 

85.0 

CCD 

94.5 

-8.4 

86.1 

SD-CI +  QC* 

87.4 

-8.3 

79.1 

Experiment:  Ref.  3 

86. Oil. 0 

Experiment:  Ref.  27 

86.7*1.5 

‘Units:  kcalmoT1. 

“Reference  4.  QC  implies  that  Davidson's  formula 
was  used  to  estimate  the  effect  of  quadruple  excita¬ 
tions. 

kcal  mol'1, 33  respectively.  The  electronic  energy  of 
the  hydrogen  atom,  using  our  basis  set,  is  -  0. 49778 
hartree.  Goddard  and  Schaefer  computed  the  zero- 
point  energies  using  vibrational  frequencies  derived 
from  force  constants  predicted  by  ab  initio  SCF  calcu¬ 
lations.9  Although  the  magnitude  of  individual  frequen¬ 
cies  differed  from  the  experimental  values,  the  zero- 
point  energy  difference  determined  from  SCF  theory 
differs  little  from  the  experimental  value. 

A  series  of  computed  dissociation  energies  is  pre¬ 
sented  in  Table  IV.  Our  self-consistent  field  calcula¬ 
tions  predict  a  dissociation  energy  (68.6  kcal  mol"1)  that 
is  much  lower  than  the  experimental  values5,34  listed  in 
the  table.  The  two  MBPT  results  (D-MBPT(4),  86.4 
kcal  mol'1);  SDQ-MBPT(4),  85. 0  kcal  mol"1)  and  the 
CCD  result  (86. 1  kcal  mol'1)  agree  well  with  the  ex¬ 
perimental  values,3,34  while  the  GS-CI  calculations  pre¬ 
dict  a  dissociation  energy  of  approximately  79  kcal 
xmol"1.5  The  present  MBPT/CCD  calculations  employ 
UHF  reference  functions  for  the  radical  products,  while 
the  SD-CI  calculations5  use  a  two-configuration  RHF- 
SCF  root  function.  In  the  latter  calculation,  an  esti¬ 
mate  of  higher  excitation  effects  is  obtained  by  using  a 
two-configuration  modification  of  Davidson’s  formula, 13 
although  the  justification  for  using  such  a  formula  with 
a  two-configuration  reference  is  not  clear. 35  For  for¬ 
maldehyde  the  net  effect  is  less  than  1  kcal  mol"1  in  any 
event.  The  linked  diagram  and  SD-CI  results  together 
demonstrate  that  high  quality  theoretical  calculations 
can  predict  energy  differences  with  good  accuracy. 

B.  CHjO-H,+CO 

Since  MBPT/CCD  results  are  size-extensive,  dis¬ 
sociation  energies  may  be  obtained  without  requiring 


TABLE  V.  Summary  of  electronic 
structure  results  for  H2  +  CO.* 


Level  of  calculation 

Total  energy 

SCF 

-113.8977 

D-MBPT<4) 

-  114.2762 

SDQ-MBPT14) 

-114.2780 

CCD 

-114.2709 

‘Units:  hartree. 


TABLE  VI.  Predicted  and  experimental  values 
for  the  dissociation  of  formaldehyde  to  molecu¬ 
lar  products.* 


Level  of  calculation 

A£j, 

SCF 

0.1 

-7.1 

-7.0 

D-MBPT(4) 

4.4 

-7.1 

-2.7 

SDQ— MBPT(4) 

3.1 

-7.1 

-4.0 

CCD 

3.9 

-7.1 

-3.2 

Experiment:  Ref. 

36 

-1.9 

•Units:  kcal  mol'1. 


supermolecule  calculations.  Table  V  contains  a  sum¬ 
mary  of  electronic  energies  for  the  supermolecule 
Hj  •  •  •  CO  obtained  by  addition.  The  zero-point  energy 
change,  computed  with  experimental  frequencies  is 
-7.1  kcal  mol'1.  Predicted  heats  of  reaction,  for 
f  =  0K,  are  gathered  in  Table  VI.  The  experimental 
value  (-1.9  kcal  mol'1)  corresponds  to  the  value  derived 
from  heat-of-formation  data  suggested  by  Benson, 33 
corrected  to  zero  degrees.  The  SCF  results  yield  a 
heat  of  reaction  that  is  more  negative  than  that  obtained 
from  either  of  the  MBPT  or  the  CCD  calculations.  The 
MBPT  and  CCD  calculations  predict  a  heat  of  reaction 
that  is  less  than  that  estimated  from  the  thermochemi¬ 
cal  data. 

In  both  dissociation  energy  estimates  described  here, 
the  high-quality  (DZP)  ab  initio  electronic  structure  cal¬ 
culations  predict  dissociation  energies  that  are  below 
the  experimentally  obtained  values.  Since  the  experimen¬ 
tal  studies  for  these  two  reactions  have  good  precision,  it 
seems  relevant  to  inquire  whether  theory  will  generally 
underestimate  dissociation  energies.  One  aspect  of  the 
computations  suggests  that  theory  will  generally  under¬ 
estimate  the  change  in  electronic  energy  in  going  from 
reactant  to  dissociation  products.  The  reason  for  this 
is  that  the  finite  basis  sets  used  in  the  molecular  calcu¬ 
lations  should  be  better  suited  to  computations  on  the 
smaller  product  fragments  than  on  the  reactant  mole¬ 
cule.  37  In  a  future  article  describing  a  series  of  hy¬ 
drogen-dissociation  reactions,  we  will  report  results 
that  suggest  that  MBPT  calculations  of  the  quality  used 
here  predict  dissociation  energies  that  are  generally 
low  by  about  1-3  kcal  mol'1.33  In  each  case  reported, 
however,  the  zeroth-order  reference  function  was  well- 
represented  by  a  single  determinant. 

V.  MOLECULAR  PRODUCTS  TRANSITION  STATES 

The  dissociation  of  formaldehyde  to  form  molecular 
products  may  occur  by  one  of  several  proposed  com¬ 
peting  pathways.  Two  possible  reaction  channels  were 
studied  by  Goddard  and  Schaefer,  leading  to  the  identi¬ 
fication  of  two  reaction  intermediates.5  One  of  these 
intermediates,  represented  schematically  in  Fig.  1(b), 
occurs  on  the  S0  potential  surface  in  the  channel  leading 
to  the  direct  formation  of  molecular  products.  The 
other  intermediate,  represented  schematically  in  Fig. 
1(c),  occurs  on  the  S0  potential  surface  in  the  channel 
leading  to  formation  of  the  metastable  species,  hy- 
droxycarbene.  Because  the  dissociation  energies  pre- 
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dieted  by  the  fourth-order  MBPT  and  the  CCD  calcula¬ 
tions  differed  by  several  kcal  mol'1  from  the  dissocia¬ 
tion  energy  predicted  by  the  SD-CI  calculations,1  it  is 
important  to  determine  whether  the  same  difference 
occurs  for  energy  barrier  predictions. 

The  calculations  reported  in  Tables  VII  and  VIII  re¬ 
fer  to  structural  parameters  determined  at  the  SCF  lev¬ 
el  using  the  gradscf  computer  codes.*4  As  noted  by 
GYS,7U>  the  results  of  the  barrier  predictions  should 
not  be  greatly  affected  by  use  of  structural  parameters 
determined  from  SCF  results,  since  the  transition 
states  are  well  described  by  a  single  determinant  wave 
function.  The  vibrational  frequencies  for  the  three 
formaldehyde  structures,  determined  using  our  (431/31) 
basis  set,  are  presented  in  Table  IX,  along  with  exper¬ 
imental  values39  and  several  other  theoretical  esti¬ 
mates.  7<b,'4°  All  calculations  were  done  using  basis 
sets  that  included  polarization  functions.  Although  the 
vibrational  frequencies  predicted  using  correlated  wave 
function  calculations  agree  better  with  experimental 
harmonic  frequencies  than  do  the  SCF-level  estimates, 
none  of  the  theoretical  methods  succeeds  in  obtaining 
good  agreement  with  experimental  values.  Yamaguchi 
et  a/. 7<b)  suggest  that  the  effects  of  higher  excitations, 
perhaps  unlinked  cluster  contributions,  have  a  signifi¬ 
cant  effect  on  the  theoretical  computation  of  frequencies. 
A  study  by  Pople  et  at.*1  supports  that  conclusion.  Cal¬ 
culations  for  12  vibrational  frequencies  using  second- 
order  Moller-Plesset  theory,  with  a  molecular  basis 
that  includes  polarization  functions  on  heavy  atoms, 
yield  results  that  differ,  on  the  average,  from  experi¬ 
mental  harmonic  frequencies  by  47.7  cm'1,  or  2.  5%. 41 

It  is  more  germane  to  the  study  of  the  barriers  on  the 
S0  surface,  however,  to  note  the  variation  in  zero-point 
energy  differences  occasioned  by  the  different  levels  of 
theory.  The  results  reported  by  GYS71*’  clearly  dem¬ 
onstrate  that  the  zero-point  energy  differences  shift 
little  when  different  levels  of  theory  are  used.  For  the 
two  barriers  reported  in  that  study,  the  SCF  and  Cl 
calculations  based  upon  the  DZ  +  P  basis  lead  to  vibra¬ 
tional  energy  correction  differences  of  0.4  and  0.  2 


TABLE  VII.  Electronic  structure  calcu¬ 
lations  for  the  molecular  products  transi¬ 
tion  state.* 


Theoretical  method  Electronic  energy 
SCF 


This  work 
DZ  ►  P“ 

DZ  *  P  Cl6 
D-MBPT(4)“ 

D-MBPT(4)d 
SIXJ— MBPT(4)d 
CCD* 

CCD* 


-113.7319 

-113.7265 

-  114.0391 
-114.1321 

-114.1310 

-114.1314 

-114.1241 

-114.1230 


•Unite:  hartree. 

“Reference  7(a). 

“Structural  parameters,  Ref.  5. 
Structural  parameters,  this  work. 


TABLE  VIH.  Electronic  structure  calcu¬ 
lations  for  the  rearrangement  transition 
state.* 


Theoretical  method  Electronic  energy 
SCF 


This  work 
DZ  +  P“ 

DZ  +  P  Cl" 

D-MBPT(4)“ 

D-MBPT(4)d 

SDQ-MBPT(4)d 

CCD“ 

CCDd 


-113.7397 

-113.7366 

-114.0481 

-114.1365 

-114.1365 

-114.1380 

-114.1293 

-114.1295 


•Units:  hartree. 

“Reference  7(a). 

“Structural  parameters,  Ref.  5. 
“Structural  parameters,  this  work. 


kcal  mol"1.  Larger  differences  occur  between  the  SCF 
and  Cl  for  the  calculations  that  did  not  include  polariza¬ 
tion  functions.  Since  our  basis  set  is  of  the  same  qual¬ 
ity  as  the  DZ  +  P  basis  used  by  GYS71*’  the  zero-point 
energy  difference  predicted  by  our  SCF  calculations 
should  be  as  accurate  as  those  obtained  by  those  authors. 
The  zero-point  energy  differences  were  computed  using 
the  frequencies  presented  in  Table  IX.  The  vibrational 
correction  for  the  molecular  products  dissociation  is 
-5.7  kcal  mol'1,  and  that  for  the  rearrangement  reac¬ 
tion  is  -4.3  kcal  mol'1.  These  values  agree  well  with 
those  determined  by  GYS71*'  and  Harding  et  at. , 30  ir¬ 
respective  of  the  theoretical  method  used  to  estimate 
the  frequencies. 

Table  X  contains  a  summary  of  the  predicted  energy 
barriers  for  the  formaldehyde  reactions.  We  report 
the  results  of  SCF,  fourth-order  MBPT  calculations  in¬ 
cluding  double-excitation  diagrams,  and  including  all 
single-,  double-  and  quadruple-excitation  diagrams  that 
contribute  at  that  order,  and  CCD  calculations.  These 
results  may  be  compared  to  the  Cl  results, 7<*’  the 
Moller-Plesset  results  of  Harding  et  al. , 39  and  to  the 
predicted  dissociation  energy  for  the  formation  of  radi¬ 
cal  products.  In  each  case,  the  SCF  calculations  pre¬ 
dict  barriers  that  exceed  those  predicted  by  the  corre¬ 
lation  methods.  Further,  the  correlated  calculations 
agree  well  with  one  another  in  the  prediction  of  the  bar¬ 
rier  heights. 

Isomerization  of  formaldehyde  to  form  hydroxycarbene 
exemplifies  the  1,  2-hydrogen  shift,  an  important  class 
of  chemical  reactions.  A  recent  review  by  Schaefer 
elucidates  the  features  of  this  reaction  and  provides  a 
review  of  the  experimental  literature. 42  The  energy 
barriers  predicted  by  the  various  correlated  methods 
range  from  85.  7  (SDQ-MBPT)  to  88.  2  kcal  mol'1 
(SD-CI). 7<*’  Each  theoretical  technique  predicts  that 
the  barrier  to  formation  of  molecular  products  exceeds 
the  rearrangement  barrier.  However,  the  correlated 
predictions  for  this  barrier  range  from  88.0  (RMP-4)30 
to  92.8  kcal  mol'1  (SD-CI), 7<t>  with  the  difference  in  the 
barrier  heights  equal  to  approximately  3.0  kcal  mol'1. 
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TABLE  IX,  Vibrational  frequencies:  CH20(S„). 


Equilibrium  structure  mode 

vi  (a,) 

v2  (a,) 

Vj  («*j) 

(ft)) 

Vf  (ft2) 

ve  (ft2) 

Experiment* 

2944 

1764 

1563 

1191 

3009 

1287 

Experiment” 

2782 

1746 

1500 

1167 

2843 

1247 

SCF 

This  work 

3191 

2004 

1659 

1343 

3113 

1368 

Ref.  7(b) 

3149 

2006 

1656 

1335 

3226 

1367 

Ref.  40 

2971.2 

1882.6 

1600.1 

1301.4 

3070.6 

1349.7 

DZtPCl  [Ref.  7(b)  1 

3074 

1869 

1596 

1243 

3155 

1306 

CEPA 

2925.6 

1797.2 

1580.9 

1220.9 

3032.3 

1308.9 

Rearrangement  transition  state 

V)  (o') 

v2  (a') 

v3  (o') 

v4  (o') 

V5  (o') 

v>6  («") 

SCF 

This  work 

3155 

2806 

1637 

1408 

264  5t 

637 

Ref.  7(b) 

3180 

2832 

1639 

1408 

2650t 

625 

DZ  *  P  Cl 

Ref.  7(b) 

3057 

2748 

1512 

1356 

2356 i 

~600d 

Molecular  products  transition  state 

Vl  (a') 

v2  (o') 

v3  (o') 

v<  (a') 

v6  (a')c 

v6  (a") 

SCF 

This  work 

3250 

2101 

1484 

819 

2295i 

1043 

Ref.  7(b) 

3243 

2092 

1526 

829 

2305) 

1024 

DZ  *  PCI 

Ref.  7(b) 

3263 

1939 

1555 

876 

2124i 

~  950d 

“Reference  39(a). 

'Reaction  coordinate. 

“Reference  39(b). 

Estimated. 

The  MBPT/CCD  results  suggest  that  the  energy  re¬ 
quirements  for  all  three  unimolecular  processes  are 
approximately  equal.  Since  it  is  often  assumed  that 
even  good  correlated  wave  function  calculations  over¬ 
estimate  barrier  heights  and  underestimate  dissociation 
energies,  the  MBPT/CCD  results  imply  that  the  rear¬ 
rangement  to  hydroxycarbene  is  the  energetically  fav¬ 
ored  reaction  pathway.  Such  a  result  agrees  with  the 
conclusions  of  Goddard  et  at. 5’7,*)  However,  basis  set 
changes  and  the  inclusion  of  triple  excitation  diagrams 
may  change  this  result. 

It  is  interesting  to  note  the  substantial  difference  be¬ 
tween  the  energy  difference  obtained  from  correlated 
calculations  and  the  results  obtained  using  just  SCF  cal¬ 
culations.  As  stated,  the  results  obtained  with  MBPT, 
CCD,  and  SD-CI  imply  that  the  energy  requirements 
for  all  three  reaction  pathways  are  almost  equal.  The 
SCF  calculations  predict  that  the  barriers  for  dissocia¬ 
tion  to  molecular  products  (95. 7  kcal  mol’1)  and  the  bar¬ 
rier  for  rearrangement  to  trans -hydroxycarbene  (99.  2 
kcal  mol'1)  differ  by  a  relatively  small  amount  of  ener¬ 
gy.  However,  the  SCF  results  imply  that  each  of  these 
barriers  is  much  larger  than  the  barrier  for  dissocia¬ 
tion  to  radical  products  (88.6  kcal  mol*1).  These  ener¬ 
gy  differences  demonstrate  the  importance  of  including 


TABLE  X.  Predicted  energy  barriers  for  formaldehyde  reac¬ 
tions.* 


Reaction 


Molecular  Radical 

Theoretical  method  1, 2-shift  products  products 


SCF 


This  work 

Ref.  7(a) 

95.7 

95.2 

99.2 

100.2 

68.6 
.  ,,b 

D-MBPT(4) 

SDQ— MBPT(4) 

86.9 

85.7 

89.0 

88.4 

86.4 

85.0 

CCD 

87.6 

90.3 

86.1 

Cl  [Ref.  7(a)) 

SD-CI 

SD-CI +QCd 

88.2 

85.9 

92.8 

88.9 

80.0 

79.1 

RMP4 

Ref.  30 

86.2 

88.0 

‘Units:  kcal  mol'*. 

”No  comparable  SCF  results  reported. 

"Includes  correction  for  quadruple  excitations  based  upon 

Davidson's  formula. 
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correlation  energy  effects  to  obtain  chemically  accurate 
results. 

An  additional  result  of  the  configuration  interaction 
study  was  that  each  decomposition  reaction  possesses  an 
energy  barrier  relative  to  the  St  band  origin.  To  de¬ 
termine  whether  the  linked-diagram  related  methods 
will  also  predict  the  existence  of  this  barrier,  we  must 
consider  the  first  excited  state  of  formaldehyde. 

VI.  THE  a 3 A"  STATE  OF  FORMALDEHYDE 

The  ultraviolet  absorption  spectrum  of  formaldehyde 
vapor  consists  of  two  band  systems:  a  strong  band  with 
apparent  origin  at  3530  A  (3. 454  eV)  and  a  weak  band  with 
apparent  origin  at  3960  A  (3. 124eV).  Walsh  predicted 
that  these  states  would  have  nonplanar  structure, 43  and 
this  was  verified  for  both  the  singlet44'45  and  the  trip¬ 
let48'47  states.  Jones  and  Coon  derived  experimental 
structures  for  the  excited  states,  including  a  specifica¬ 
tion  of  the  out-of-plane  bending  angle  by  a  distribution 
function  for  each  vibrational  level  of  the  bending  mode.48 
Recent  theoretical  studies  of  the  excited  states  include 
ab  initio  calculations  by  Garrison,  Schaefer,  and  Les¬ 
ter,  10  Bell,8  and  Goddard  and  Schaefer, s  as  well  as  equa- 
tions-of-motion  calculations  of  the  vertical  excitation 
energies. 17,18  This  report  complements  and  extends 
Bell’s  work. 

Recent  calculations5'8  used  excited  state  structural 
parameters  determined  by  optimizing  the  electronic 
energy  at  the  SCF  level.  This  study  reports  structural 
parameters  for  the  a3 A”  state  determined  by  optimizing 
the  energy  at  the  correlated  level.  Fourth  order,  dou¬ 
ble-excitation  MB PT  calculations,  with  core  orbitals 
frozen,  were  performed  for  structures  obtained  by  vary¬ 
ing  the  CH  and  CO  bond  lengths,  the  HCH  angle,  and  the 
out-of-plane  bending  angle.  The  CH  bond  length  and  the 
HCH  angle  were  reoptimized  after  determination  of  the 
optimal  out-of-plane  bending  angle.  We  report  the  op¬ 
timal  structural  parameters  in  Table  XI. 

The  carbon-hydrogen  bond  length  (1.085  A)  reported 
here  is  slightly  less  than  the  experimental  value 
(1.0962  A)  while  the  carbon-oxygen  bond  length  (1.327 
A)  exceeds  the  experimental  value  (1.  307  A).  The  op¬ 
timum  HCH  and  out-of-plane  bending  angles  (118°  and 
37°  12')  agree  well  with  the  experimental  values  (118° 
and  37° 54').  Comparison  of  these  theoretical  results 
with  Bell’s  structure8  shows  that  the  inclusion  of  cor- 


TABLE  XII.  Adiabatic  excitation  energies, 
S0  ~Tt. 


AE  (eV) 

A E  (kcal  mol*1) 

Experiment* 

3.124 

71.9 

D-MBPT(4) 

3.1455 

72.5 

SDQ-MBPT14) 

3.1436 

72.5 

CCD 

3.093 

71.3 

SD-CI“ 

3.01 

69.5 

SD-CI“ 

2.664 

61.4 

scf“ 

2.183 

50.3 

SCF,  this  work 

1.912 

44.1 

“Reference  32. 
“Reference  5. 
“Reference  6. 
“Reference  10. 


relation  effects  improves  the  agreement  between  theory 
and  experiment. 

We  report,  in  Table  XII,  three  theoretical  results  for 
the  adiabatic  excitation  energy  corresponding  to  S0-  T,. 
Both  the  D-MBPT(4)  and  SDQ-MBPT(4)  calculations 
yield  72.  5  kcal  mol"1  for  the  excitation  energy,  while 
the  CCD  calculations  yielded  a  value  of  71.3  kcal  mol"1. 
These  results  may  be  compared  to  the  near  Hartree- 
Fock  value,  50.3  kcal  mol'1, 10  two  SD-CI  values  67.4 
kcal  mol*1 5  and  68.9  kcal  mol"1, 49  and  a  value  obtained 
by  including  the  unlinked  cluster  correction  of  quadruple- 
excitation  effects  69.  5  kcal  mol"1. 5  All  the  theoretical 
values  compare  well  with  the  experimental  value  71.9 
kcal  mol"1. 32  None  of  the  theoretical  values,  however, 
include  any  correction  for  zero-point  energy  differences 
that  may  contribute  to  the  observed  value. 

Although  no  effort  has  been  made  to  determine  the 
structure  of  formaldehyde  in  the  A  ‘A"(S,)  state,  ver¬ 
tical  excitation  energies  are  computed  for  both  the 
S0-Si  andS0—  Tj  transitions.  Since  the  UHF  wave 
function  for  the  S,  state  was  contaminated  by  higher- 
order  spin  states,  particularly  T,,  the  computed  exci¬ 
tation  energy  for  this  state  is  expected  to  be  less  ac¬ 
curate  than  that  for  the  T,  state.  Nonetheless,  vertical 
excitation  energies  for  both  transitions  are  presented 
in  Table  Mil,  along  with  other  theoretical  predictions 
and  experimental  estimates.  The  excitation  energy  for 
S0-  T,  is  less  than  that  for  S0-S„  although  the  latter 
is  probably  lower  than  would  be  obtained  were  the  ref- 


TABLE  XI.  Formaldehyde  structural  parameters,  a3 A" 
state. 


rcn  * 

rCO* 

HCH 

<t> 

Experiment8 

1.0962 

1.3069 

118” 

37”  54' 

Theory 

DZ  P-SC  F“ 
D-MBPT(4)-FC 

1.0779 

1.0848 

1.3404 

1.3271 

119.67” 

118” 

38”  30' 

37”  12' 

“Units:  angstroms. 
“Reference  48. 
“Reference  8. 


TABLE  XIII.  Vertical  excitation 
energies  for  formaldehyde.* 


s0-r, 

S0  “Si 

D-MBPTI4) 

84.2 

89.  J 

MRD-CI“ 

78.6 

87.9 

EOM“ 

79.8 

93.2 

Experiment4 

83.0 

98.2 

“Units:  koalmol*1. 

“Reference  50. 

“Reference  18. 

“Estimated  excitation  energies, 
Ref.  51. 
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TABLE  XIV.  Energy  prediction  comparisons  for  photochemical  processes.* 


Reference  5 

Reference  7(a) 

SDQ-MBPTI4) 

CCD 

Experiment 

s0-r, 

69.5 

72.5 

71.3 

71.9* 

O. 

~  79. 

85.0 

86.1 

86. 0±16 

£ o  <H2CO) 

87.5 

88.9 

88.4 

90.3 

£» 

83.9 

85.9 

85.7 

87.6 

s„— s, 

78(est.) 

81. 0(est. ) 

79. 8(est. ) 

80.5* 

‘Reference  32. 

‘Reference  3. 

“Reference  1. 

erence  function  not  spin-contaminated.  The  difference  in 
the  excitation  energy  predictions  is  much  less  than  the 
difference  obtained  from  the  equation  of  motion  calcu¬ 
lations17  or  estimated  from  experimental  data. 51  How¬ 
ever,  it  should  be  noted  that  MRD-CI  calculations  pre¬ 
dict  a  significantly  smaller  excitation  energy  difference 
(9. 3  kcal  mol'1)50  than  is  estimated  by  Chutjian  (15.  2 
kcalmol'1).51  Since  the  difference  between  the  two  adia¬ 
batic  excitations  is  8. 1  kcal  mol'1,  it  seems  best  to  con¬ 
clude  that  the  S0-St  vertical  excitation  energy  remains 
uncertain. 

VII.  DISCUSSION 

The  present  calculations  demonstrate  that  linked-dia- 
gram  related  computational  methods  predict  relative 
energies  of  small  molecular  systems  to  high  accuracy. 
For  the  results  that  can  be  compared  with  experiment, 
the  dissociation  energy  for  formation  of  radical  prod¬ 
ucts  and  the  Sa  -  Tx  adiabatic  excitation  energy,  the 
SDQ-MBPT(4)  and  the  CCD  results  agree  well  with  the 
experimental  values. 

It  is  instructive  to  review  results  pertinent  to  photo¬ 
chemical  dissociation  and  rearrangement  of  formalde¬ 
hyde.  We  collect,  in  Table  XIV,  the  excitation  energy, 
dissociation  energy  (D„),  and  reaction  barrier  (£0)  re¬ 
sults  obtained  in  this  study  and  in  work  by  Goddard  and 
co-workers.  5,7<*)  Note  that  the  values  for  S0-S,  are 
estimates  based  on  the  theoretical  results  for  the  S0~  T, 
excitation  energy,  and  the  observed  T1-S1  separation. 

In  analyzing  the  Cl  results,  the  authors  noted  that  high 
level  theoretical  calculations  underestimate  dissocia¬ 
tion  energies  and  slightly  overestimate  rearrangement 
barriers.52’53  If  it  is  assumed  that  the  Cl  calculations 
underestimate  the  dissociation  energy  by  3  to  5  kcal 
x  mol*1,  and  overestimate  the  energy  barriers  by  a 
similar  amount,  then  the  Cl  results  imply  that  the  en¬ 
ergy  requirements  for  the  radical  products  and  molecu¬ 
lar  products  dissociation  pathways  are  similar.5  Fur¬ 
thermore,  the  energy  requirements  suggest  that  the  hy- 
droxycarbene  intermediate  cannot  be  ruled  out.  The 
subsequent  analytic  gradient  Cl  calculations  support 
these  results,  but  each  of  the  energy  barriers  increases 
slightly. 7U>  Finally,  those  authors  noted  that  the  bar¬ 
riers  exceed  the  estimated  adiabatic  excitation  energy 
for  S„-S,.  Thus,  there  all  small  barriers  above  the 
S,  origin  to  all  three  photochemical  processes. 

Examinations  of  the  SDQ-MBPT/CCD  results  reveals 
no  significant  difference  from  the  results  of  the  Cl  cal¬ 
culations,  even  though  our  basis  set  is  somewhat  larger 
than  the  GS  DZ  +  P  basis.  However,  our  SCF  result  for 


the  molecular  products  transition  state  shows  an  energy 
lowering  of  3.4  kcal/mol  compared  to  the  DZ  +  P  re¬ 
sults,  suggesting  that  basis  set  changes  may  still  be 
significant  in  this  problem.  We  cannot  make  realistic 
estimates  of  the  errors  associated  with  the  barrier  cal¬ 
culations,  since  we  have  no  data  for  comparison.  We 
assume  that  these  calculations  overestimate  the  energy 
barriers  by  2  to  4  kcalmol'1.  With  these  assumptions, 
the  results  of  these  calculations  lead  to  the  same  conclu¬ 
sions  obtained  by  GS5  and  Goddard  et  al. ,u>  However, 
recent  work  by  Frisch  et  al.**  suggests  that  basis  set 
effects  and  triple  excitations  lower  the  computed  bar¬ 
riers  by  "8  kcal/mol,  which  is  below  the  origin  of  the 
H2CO  excited  singlet  state. 
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